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The optimized perturbation theory (OPT) at finite tem- 
perature (T) recently developed by the present authors is re- 
viewed by using O(N) 4> 4 theory with spontaneous symmetry 
breaking. The method resums automatically higher loops (in- 
cluding the hard thermal loops) at high T and simultaneously 
cures the problem of tachyonic poles at relatively low T. We 
prove that (i) the renormalization of the ultra-violet diver- 
gences can be carried out systematically in any given order of 
OPT, and (ii) the Nambu-Goldstone theorem is satisfied for 
arbitrary N and for any given order of OPT. 



I. INTRODUCTION 

Naive perturbation theory is known to break down at 
finite temperature (T) . The two reasons are the existence 
of hard thermal loops (HTL) at high T @ and the emer- 
gence of tachyonic poles at relatively low T ||. If one 
adopts self-consistent resummation methods, HTL can 
be summed and tachyonic poles can be removed. How- 
ever, most of the self-consistent methods proposed so far 
have difficulties of renormalization H and/or the viola- 
tion of the Nambu-Goldstone theorem [Q at finite T. 

In this talk, we show that a new loop-wise expansion 
at finite T recently developed by the present authors || 
can solve these problems. 

Our starting point is the optimized perturbation the- 
ory (OPT) which is a generalization of the mean-field 
method jq| and is known to work in various quantum 
systems M . Its application to field theory at finite T has 
been first considered by Okopihska || and Banerjee and 
Mallik We further develop the idea and prove the 
renormalizability and the Nambu-Goldstone (NG) the- 
orem in O(N) </> 4 theory at finite T order by order in 
OPT. 

The organization of this talk is as follows. In section 
II, we introduce a loop-wise expansion on the basis of 
OPT. The renormalization of UV divergences and the 
realization of the NG theorem in this method are also 
discussed. In section III, OPT is applied for the 0(4) 



a model which is a low energy effective theory of QCD. 
Summary and concluding remarks are given in Section 
IV. 



II. OPTIMIZED PERTURBATION ATT/O 
A. Hard thermal loops and tachyonic poles 

Let us illustrate, by using </> 4 theory, the reason why 
the naive perturbation theory at finite T breaks down; 
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We first consider the case fj, 2 > 0. The lowest order self- 
energy diagram Fig.0 (A) is <3(AT 2 ) at high T. However, 
Fig.[l] (B) is 0(XT 2 x — ). Furthermore, higher powers 
of T/fi arise in higher loops; e.g. the n-loop diagram in 
Fig.| (C) is 0(\ n T 2n ~ 1 / f i 2n - 3 ). Thus, the validity of the 
perturbation theory breaks down when T > fi/X because 
the higher order diagrams are larger than lower ones. 
Therefore, one should at least resum cactus diagrams to 
get sensible results at high T Physics behind this 
resummation is the well-known Debye screening mass in 
the hot plasma. 

When /i 2 < and the system has spontaneous symme- 
try breaking (SSB), the naive perturbation shows another 
problem. The tree-level mass tuq in this case is defined 
as 



2 ' 



(2) 



where £(T) is the thermal expectation value of <p. As 
T increases, £ decreases. Then m§ becomes negative 
(tachyonic) even before the critical temperature T c is 
reached. If this happens, the naive perturbation using 
the tree-level propagator does not make sense and cer- 
tain resummation should be carried out [|l0|. Note that, 
for T < T c , there is no reason to believe that only the 
cactus diagrams shown in Fig.l are dominant; there ex- 
ists a three-point vertex A£<^> 3 which is not negligible for 
T~£(T). 



*Talk presented at "Thermal Field Theories and Their Ap- 
plications", (Regensburg, Germany, August 10-14, 1998) 
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B. Problems in self-consistent resummation methods 

Self-consistent resummation method is a procedure to 
improve perturbation theory at finite T and to avoid the 
problems in Sec. 



II A 



However, the method has other 

difficulties §,[§. 

In the naive perturbation theory, there arises no new 
UV divergences at T ^ because of the natural cutoff 
from the Boltzmann distribution function. Therefore, all 
the UV divergences at finite T are canceled by the counter 
terms prepared at T — jnj. 

On the other hand, in self-consistent methods at T ^ 
0, the situation is not that simple: In fact, the tree- 
level propagators have T-dependent mass (such as m(T) 
in the above) which contains higher loop contributions 
through the self-consistent gap-equation ||]. This leads 
to a necessity of T-dependent counter terms which are 
sometimes introduced in ad hoc ways. 

Another problem is the violation of the Nambu- 
Goldstone (NG) theorem: In many of the self-consistent 
methods, resummation with keeping symmetry is a non- 
trivial issue, and the NG theorem is often violated. 



C. New resummation method 

For theories with SSB, loop-expansion rather than the 
weak-coupling expansion is relevant, since one needs to 
treat the thermal effective potential. Therefore, we de- 
veloped an improved loop-expansion at finite T for the 
purpose of resummation |^] . The method keeps the renor- 
malizability and guarantees the Nambe-Goldstone theo- 
rem order by order at finite T . 

In the following, we divide our resummation procedure 
into three steps and apply it to <j) A theory. The case for 
O(N) (f) 4 theory will be discussed in Sec. HE. 



We start with the thermal effective action with an ex- 
pansion parameter "<5" : 



In 



|exp 



l/T 



d n x [£( 



where J = — dT[ip]/d(p and J^ T d n x : 
If we explicitly write ft in eq. 

rh/T 



f^drfd^x. 
it appears as 

j- Jq" /J d A xC. Therefore, the loop-expansion by 5 at finite 
T does not coincide with the ft, expansion. The expansion 
by 8 should be regarded as a steepest descent evaluation 
of the functional integral. 

Step 1 

Start with a renormalized Lagrangian with counter 
terms 



J<j 



i 



A 

4!< 



-Bfi 



C^ + D f i 4 . (4) 



Here we have explicitly written the argument /z 2 in £ for 
later use. The MS scheme with the dimensional regular- 
ization is assumed in (Q). Just for notational simplicity, 
the factor k^ 4- ™) to be multiplied to A is omitted (k is 
the renormalization point). 

The c- number counter term -D/i 4 , which was not con- 
sidered in [||, is necessary to make the thermal effective 
potential finite. Also, it plays a crucial role for renormal- 
ization in OPT as will be shown in Sec. 



II D 



The renormalization constants are completely fixed at 
T = and fi 2 > 0, in which A, B, C and D are expanded 
as 



(5) 



Note that (i) the coefficients (a/, 6;, q, di) are independent 
of /i 2 , since we use the mass independent renormalization 
scheme, (ii) the UV divergences in the symmetry broken 
phase (/i 2 < 0) can be removed by the same counter 
terms determined for /i 2 > [[l2 13|, and (iii) A,B,C 
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and D are independent of T by definition. 

The relations of A, B, C and D with the standard 
renormalization constants are A = Z — 1, B = Z^Z — 1 
and C = Z\Z 2 — 1, where Z's are defined by 4>q = \[Z^>^ 
Ao = Z\\ and /i 2 = Z^fi 2 with suffix indicating un- 
renormalized quantities. 



Step 2 

Rewrite the Lagrangian by introducing a new mass 
parameter m 2 following the idea of OPT M: 



2 2 / 2 2\ _ 2 

fx r = m — [m — fi ) = m — X- 



(6) 



This identity should be used not only in the standard 
mass term but also in the counter terms [3, which is 
crucial to show the order by order renormalization in 
OPT: 



(3) 



£{4; ^ 2 



C x + £ c 



(7) 



A-ip0) 2 -mV]-^4 + l x0 2 (8) 
C c = \A{d<t>) 2 - l -B{m 2 - xW ~ jf<f 



+D(m 2 - x? 



(9) 



A, B, C and D in C c were already determined in Step 1. 

On the basis of eq.(0), we define a "modified" loop- 
expansion in which the tree-level propagator has a mass 
m 2 + \ip 2 /2 instead of /j, 2 + X(f 2 /2. Major difference be- 
tween this expansion and the naive one is the following 
assignment, 



to 2 = 0(S°), X = 0(6). 



(10) 



The physical reason behind this assignment is the fact 
that x reflects the effect of interactions. If one makes an 
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assignment, m 2 = 0(6°), X — 0(6°), the modified loop- 
expansion immediately reduces to the naive one. 

Since eq.(Q) is simply a reorganization of the La- 
grangian, any Green's functions (or its generating func- 
tional) calculated in the modified loop-expansion should 
not depend on the arbitrary mass to if they are calculated 
in all orders. However, one needs to truncate perturba- 
tion series at certain order in practice. This inevitably 
introduces explicit to dependence in Green's functions. 
Procedures to determine to are given in Step 3 below. [] 

To find the ground state of the system, one should look 
for the stationary point of the thermal effective potential 
defined by V(p 2 ) = — T[tp 2 = const.]/ J* d*x. As men- 
tioned above, V calculated up to L-th loops VL(f 2 ;m) 
has explicit m-dependence. Thus the stationary condi- 
tion reads 



dV L (ip 2 -m) 
dip 



= 0, 



(11) 



where derivative with respect to ip does not act on m 
by definition. Eq.([ll]) gives a stationary point of Vl for 
given to. 

Step 3 

The final step is to find an optimal value of to by im- 
posing physical conditions a la Stevenson jlj] such as the 
following. 

(a) The principle of minimal sensitivity (PMS): this 
condition requires that a chosen quantity calculated 
up to L-th loops (Ol) should be stationary by the 
variation of m: 



dO L 
dm 



0. 



(12) 



(b) The criterion of the fastest apparent convergence 
(FAC): this condition requires that the perturba- 
tive corrections in Ol should be as small as possible 
for a suitable value of m. 



Or 



L -i 



0, 



(13) 



where I is chosen in the range, 1 < I < L. 



The above conditions reduce to self-consistent gap 
equations whose solution determine the optimal param- 
eter m for given L. Thus m becomes a non-trivial func- 
tion of ip, A and T. This together with the solution 



1 One may generalize Step 2 by adding and subtracting 
ao(do4>) 2 1 ai(di(f>) 2 and 7<^) 4 with ao, ol\ and 7 being finite 
parameters to be determined by the PMS or FAC conditions 
(see Step 3). The renormalizability can be also shown to be 
maintained in this case. However, we will concentrate on the 
simplest version («o,i = 7 = 0) in the following discussions. 



of ( |1 1| ) completely determine the thermal expectation 
value £(T) = (</>) as well as the optimal parameter m(T). 
Through this self-consistent process, higher order terms 
in the naive loop expansion are resumed. 

The choice of Ol in Step 3 depends on the quantity 
one needs to improve most. To study the static nature 
of the phase transition, the thermal effective potential 
Vl[}P 2 \ m) is most relevant. Applying the PMS condition 
for Vl reads 



dV L {p> 2 ;m) 
dm 



= 0, 



(14) 



which gives a solution m — m(ip). This can be 
used to improve the effective potential at finite T as 
V L (p 2 ;m) -► V L {<p 2 ;m(tp)). Also, £(T) and m{T) are 
obtained by solving ( pd| ) together with (Q . In this case, 
the following relation holds: dV(p 2 ; m(p)) / d<p\ v —^ = 
dV{<p 2 ;m{(p))/d<p\ v= ^. 

To improve particle properties at finite T, it is more 
efficient to apply PMS or FAC conditions directly to the 
two-point functions. We will use FAC for the one-loop 
pion self-energy in Section III to show its usefulness. 



D. Renormalization in OPT 

We now prove the order by order renormalization in 
OPT. Let us first rewrite eq.(0) as 



£(0;m 2 ) = £{4>;m 2 ) 
+\x4> 2 + 



\bx4> 2 



D X - 2Dm 2 x 



(15) 



The UV divergences arising in the perturbation theory 
are classified into two classes: The divergences in the 
Green's function generated by £(0; m 2 ), and the diver- 
gences obtained by the multiple insertion of (l/2)x</!> 2 to 
the Green's function generated by £(</>; to 2 ). 

Since we use the symmetric and mass independent 
renormalization scheme (such as the MS scheme), any 
divergences in the first class are renormalized solely 
by the coefficients A, B, C and D in £(0;/i 2 ). Al- 
though T-dependent divergences appear because of the 
T-dependent "resumed" mass to 2 (T), they are prop- 
erly renormalized away since the counter terms (such 
as Bm 2 <p 2 , Dm 2 and Dm 4 ) also acquires T-dependence 
through to 2 (T). In other words, the divergences arising 
from the resumed propagator is removed by the resumed 
counter terms. (See also, footnote 1.) 

The divergences in the second class can be shown to be 
removed by the last three counter terms in (|l5|). (Note 
that B and D are already fixed in Step 1, and we do 
not have any freedom to change them.) This is obviously 
related to the renormalization of composite operators. 
In fact, the standard method jl6| tells us that necessary 
counter terms are written as 
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2 y ^ 



i)x4> 2 + A 2 x 2 + Aix- 



(16) 



Here Z^p. is the renormalization constant for the compos- 
ite operator </> 2 , and removes the divergence in Fig.|2|(A). 
A2 and Ai are necessary to remove the overall diver- 
gences in Fig.||(B) and in Fig p|(C), respectively. 

Now, one can prove that ( |l6|) coincides with the last 
three terms in CfL5|) : 



ZZ~} - 1 



B, A 2 = D, and A a = -2Dm 2 . (17) 



The first equation is obtained by the definition B 
Zf+Z — 1 and an identity 



z <p — Z n 



(18) 



The overall divergence of the vacuum diagram with no 
external-legs is removed by the c-number counter term 
Dm 4 in C((j>; m 2 ). Therefore, the last two equations in 
( |l7| ) are obtained as 



d 



2A, 



y dm 2 
(A) 2 [Dm 4 



2D. 



(19) 
(20) 



Eq. ( |l~7[ ) shows clearly that all the necessary counter terms 
in OPT are supplied solely by the original Lagrangian 
£(0;/i 2 ). Thus, we can carry out renormalization or- 
der by order even within the self-consistent method. For 
more detailed proof of the relations JT^), see Appendix 
A of ref. §. 

Three comments are in order here: 

(i) Because the renormalization is already carried out 
in Step 2, one obtains finite gap-equations from the 
beginning in Step 3. Our procedure "resummation 
after renormalization" has many advantages over 
the conventional procedure "resummation before 
renormalization" where UV divergences are hoped 
to be canceled after imposing the gap-equation. 
The difference between the two is prominent es- 
pecially in higher order calculations. 

(ii) The decomposition (|J) should be done both in the 
mass term and the counter terms. This guaran- 
tees order by order renormalization in our modified 
loop-expansion in any higher orders. (In ref. 
the renormalizability was checked up to the two- 
loop level in the 4 theory at high T.) On the 
other hand, if one keeps the original counter term 
(l/2)B[i 2 (j) 2 + Dfj, 4 without the decomposition (||), 
L-loop diagrams with L > M must be taken into 
account to remove the UV divergences in the M- 
loop order. This is an unnecessary complication 
due to the inappropriate treatment of the counter 
terms. (See e.g. ref. [O] which encounters this 
problem). 



(iii) As far as we stay in the low energy region far below 
the Landau pole, we need not address the issue of 
the triviality of the 4 theory |Q: Perturbative 
renormalization in OPT works in the same sense 
as that in the naive perturbation. 



E. Nambu-Goldstone theorem in OPT 

The procedure and the renormalization in OPT dis- 
cussed above do not receive modifications even if the La- 
grangian has global symmetry. For O(N) <fi 4 theory, one 
needs to replace 4> and 4> 2 by <j> = ($0, <pi, ■ ■ •, 4>n-i) and 
4> 2 respectively in all the previous formulas. 

In the symmetry broken phase of such theory, the 
Nambu-Goldstone (NG) theorem and massless NG 
bosons are guaranteed in each order of the modified loop- 
expansion in OPT for arbitrary N . To show this, it is 
most convenient to start with the thermal effective po- 
tential V((f 2 ). By the definition of the effective potential, 
V(ip 2 ) has manifest O(N) invariance if it is calculated in 
all orders. 

In OPT, V calculated up to L-th loops Vl(0 2 ;iti) 
has also manifest O(N) invariance, because our decom- 
position (^|) used in ([?]) does not break O(N) invari- 
ance. Once Vl has invariance under the O(N) rotation 
(ipi — ► (fi + i9 a T°-(fj), the immediate consequence is the 
standard identity: 



dV L (<p 2 ;m) 
dtpj 



rpa 

ji ~ 



d 2 V L (0 2 -m) 



(21) 



with T a being the generator of the O(N) symmetry. 
Eq.(pTl) is valid for arbitrary L, m and N. 

At the stationary point where the l.h.s. of (|l]) van- 
ishes, there arises massless NG bosons for T^ipk 7^ 0, 
since the r.h.s. of ( fj"l| ) is equal to — V^ 1 (0)T°- k 4>k where 
T>ij(0) is the Matsubara propagator at zero frequency 
and momentum calculated up to L-th loops. Thus the 
existence of the NG bosons is proved independent of the 
structure of the gap-equation in Step 3. 

Now, let us show an example of the unjustified approx- 
imations leading to the breakdown of the NG theorem. 
Suppose that we make a general decomposition such as 



I^Sij = ///;, 



(m 2 j 



2 r 



), 



(22) 



with mfj m 2 6ij. This leads to an O(N) non-invariant 
effective potential, and the relation (|2|) is not guaranteed 
in any finite orders of the loop-expansion. For example, 
when the O(N) symmetry is spontaneously broken down 
to 0(N — 1), one may be tempted to make a decompo- 
sition m 2 j = to 2 , (i = j = 0), mjj — to 2 (i = j ^ 0), and 
rnfj — (i 7^ j) to impose self-consistent conditions for 
the radial mode and the NG mode. However, the effec- 
tive potential does not have O(N) symmetry in this case 
and eq.(pl|) does not hold. 
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III. APPLICATION TO 0(4) a MODEL 

Let us apply OPT to the 0(4) a model. The model 
shares common symmetry and dynamics with QCD and 
has been used to study the real-time dynamics and criti- 
cal phenomena associated with the QCD chiral transition 
PPL 



A. Parameters at T = 



The 0(4) a model reads 



1 



A{d(f>Y 



1 



-Bp, 



2j2 



Y + ha 
A 



C{ft) 2 +Dp\ (23) 



2 -v r, 2 r T 4! 

with <j) = (a, 7?) . ha is an explicit symmetry breaking 
term. 

A, B, C and D in the one-loop order are 



A = Q, B 



A 1 



167T 2 £ ' 



C 



A 1 



D 



1 1 



i67r 2 r 



(24) 



where 4 = — 7 + log(47r), with 7 being the Euler 
constant. 

When SSB takes place (p 2 < 0), the replacement a — > 
(7 + £ in cq.(]23) leads to the tree- level masses of a and 7r; 



2 1 ^ ,-2 2 2 1 " ^2 
/i + ~£ , m 07r = /x + -£ 



(25) 



The expectation value £ at T = is determined by the 
stationary condition for the standard effective potential 
dV^/dcpj = 0. 

Later we will take a special FAC condition in which m 2 
deviates from p 2 only at T ^ 0, so that the naive loop- 
expansion at T = is valid. The renormalized couplings 
/j 2 , A and h can thus be determined by the renormaliza- 
tion conditions in the naive loop-expansion at zero T such 
as (i) m n = 140 MeV, (ii) f n = 93 MeV, and (iii) tt-tt 
scattering phase shift jlj| . 

Instead of (iii), one may adopt m£ eafc (the peak posi- 
tion of the spectral function in the a channel) . We take 
this simplified condition with three possible cases: m£ eafe 
=550 MeV, 750 MeV and 1000 MeV. mf ai = 550 MeV 
and 750 MeV are consistent with recent re-analyses of 
the 7T-7T scattering phase shift plj. 



m peak ( Me y) 


p 2 (MeV 2 ) 


A 


h (MeV a ) 


k (MeV) 


550 


-284 2 


73.0 


123 3 


255 


750 


-375 2 


122 


124 3 


325 


1000 


-469 2 


194 


125 3 


401 



TABLE I. Vacuum parameters corresponding to rri^ ak 
550, 750, 1000 MeV 



We still have a freedom to choose the renormaliza- 
tion point n. Instead of trying to determine optimal k 
by the renormalization group equation for the effective 
potential p2[ , we take a simple and physical condition 
7710^=771^=140 MeV. This choice has two advantages: (a) 
One-loop pion self-energy E T (/c 2 ) vanishes at the tree- 
mass; S^(fc 2 = to 2 ,^) = 'E 7r (k 2 = to 2 ) = 0, and (b) the 
spectral function in the a channel starts from a correct 
continuum threshold in the one-loop level. Resultant pa- 
rameters are summarized in TABLE I. 

In Fig^, the spectral functions p a and p^ at T = 0, 
namely the T = limit of eq.p^) d efined below, are 
shown as a function of ^fs = Vw 2 — k 2 . In the it channel, 
there are one particle pole and a continuum starting from 
the threshold ^Jsth = moir+moa- y/sth is the point where 
the channel n + a opens. In the a channel, the spectral 
function starts from the threshold 2mo„ = 280 MeV and 
shows a broad peak centered around yfs = mg eafe . The 
large width of a is due to a strong a-lix coupling in the 
linear a model. The corresponding cr-pole is located far 
from the real axis on the complex s plane. 

Here we show the definition of the spectral function at 
finite T: 



p*(w,k;T) = --ImD£(w,k;T), 
where D? is the retarded correlation function 



(26) 



Df(u,k;T) = -i I d 4 xe^^)([^,x), 0(0,0)]), (27) 
with (■) being the thermal expectation value. 

B. Application of OPT 

Now let us proceed to Step 2 in OPT and rewrite 
eq.(E3) as 



\[{d$) 2 

— Bm 2 ^ 2 
2 Y 



(28) 



Since X ( = to 2 — p 2 ) is already a one- loop order, we 
have neglected the terms proportional to B\, D\ 2 and 
D\ which are two-loop or higher orders. 

When SSB takes place (cf> — ► <f> + (p), the tree-level 
masses to be used in the modified loop-expansion read 



2 A _ 

TO + -If, 



2 A „ 
TO + -ip. 

6 



(29) 



Since to 2 will eventually be a function of T, the tree- 
masses running in the loops are not necessary tachyonic 
at finite T contrar y to the naive loop-expansion (see the 
discussion in Sec. [I A ). 
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The thermal effective potential V(0\ to) is calculated 
in the standard manner except for the extra terms pro- 
portional to x- The effective potential in the one- loop 
level reads 



V(<p;m) = ^mV 2 
1 



-i 





64tt 2 
T 



'Oct 



In 



hip 



'Oct 



s 2 e 3/2 



3m 0jr In 



'Ott 



2„3/2 



discussed in Sec. II E can be explicitly checked by com- 
paring eq.(^Tj) and the inverse pion-propagator at zero 
momentum [D%(0, 0; T)]" 1 . 

Let us impose the FAC condition on E^?. Since we 
chose a renormalization condition too. = m x = 140 MeV 
at zero T, one may be tempted to adopt the following 
condition at finite T: 



(30) 



d 3 k 



(27T) 



ln(l- e - E ^ T )+3 1n(l-e- E "/ T ) 



where = , / k 2 



Although this has the similar 

structure with the standard free energy in the naive loop- 
expansion, the coefficient of the first term in the r.h.s. of 
( |30| ) is /J 2 instead of m 2 . This is because we have extra 
mass-term proportional to \ m the one-loop level. The 
stationary point £ is obtained by 



dV(ip; to) 



d<pi 



0. 



(31) 



<?=(f,o) 



Since the derivative with respect to £ does not act on 
to, this gives a solution £ as a function of T and to. By 
imposing another condition on m (Step 3), one eventually 
determines both £ and to for given T. 



E*(w = m 0w> 0;T) = 0. 



(34) 



.However, eq.(|34|) does not guarantees that to 2 is real, 
since the l.h.s. of eq.(|34|) receives an imaginary part due 
to the Landau damping. To avoid this problem, we take 
a hybrid condition: 



E* 1 (w = too., 0) + = 0, 0; T) = 



ii/ 



(35) 



Note that the external energy is set to be zero in the In- 
dependent part.0 Because the second term in the l.h.s. 
vanishes at T = 0, the solution of eq.(|35|) at T = be- 
comes 



to 2 (T = 0) = n 2 



(36) 

does not spoil the 



Therefore, the FAC condition i 
naive-loop expansion at T = 0. 

At high T (£(T) ~ 0), the following analytic solution 
is obtained as far as T 2 3> to 2 ; 



C. FAC condition for m 

To resum the hard thermal loops, the PMS condi- 
tion for the effective potential requires 2-loop calculation, 
while the FAC condition for the self-energy requires only 
1-loop calculation. Therefore, we adopt the latter condi- 
tion here to determine to 2 . 

The retarded self-energy E^ (defined by [-D^ ] -1 = s — 
to 2 ,^ — E;?) is related to the 11-component of the 2x2 
self-energy in the real-time formalism ]24j ]; 

ReE«( W , k; T) = Re^ 1 (w, k) + S* 1 (w, k; T)} (32) 
ImE«( W , k; T) = tanh(^) lm{I%(u>, k) + Ej, 1 (w, k; T)}. 

Here Ej, 1 (w, k; T) is a part having explicit T-dependcncc 
through the Bose- Einstein distribution, while Si 1 (a;, k) is 
a part having only implicit T-dependence through m(T) 
and £(T). In the one-loop level, Ej, 1 can be calculated 
only by the 11-component of the free propagator, 



m 2 (T) = ^ 2 + ±T 2 , 



(37) 



which implies that the Debye screening mass at high T 
is properly taken into account. For realistic values of A 
in TABLE I, the condition T 2 » to 2 is not well satisfied 
and eq.(|5|) should be solved numerically. 

For intermediate values of T, eq.(|35|) can effectively 
sum not only the contributions from the diagrams in 
Fig.0(a, b, h, i), but also from those in Fig.[|(c, d, j). 
Thus, OPT can go beyond the cactus approximation 
which sums only the diagrams in Fig.^(a, b, h, i). 

Two remarks are in order here. 

(i) By eq. (|35|) , only the k-independent part of the 
self-energy is resumed. If one needs to resum k- 
dependent part too, one must introduce to 2 which 
depends both on T and k and impose k-dependent 
FAC or PMS conditions 



iDl x Ak 2 -T) 



k 2 - 



2irn B 5{k 2 - m 2 ), (33) 



"0<t> 



with n B = [e" /T - 

One-loop diagrams in OPT for SV 1 are shown in Fig. 
[I]. Their explicit forms are given in |5j . The NG theorem 



Eq. ( |35| ) can be formulated in a covariant way by introduc- 
ing the four-vector n M characterizing the heat bath. The first 
term of the equation is only a function of k 2 ^ because it is T 
independent. The second term is a function of ki and k^n^ . 



Therefore, the condition reads E^fc 
0, fc M n M = 0;T) = 0. 



m^) + E^(fc,1 
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(ii) In ref. |T^|, it has been studied the convergence 
properties of the free energy at high T with a vari- 
ational condition equivalent to the PMS condition 
here. Although the approach has a problem of 
renormalization as we have already mentioned, the 
result is suggestive in the sense that the optimized 
loop expansion has much better convergence prop- 
erties than the loop expansion based on the hard 
thermal loops [p5|. Better understanding of the 
convergence properties both in PMS and FAC con- 
ditions is an important future problem. 



D. Behavior ofm(T), m 0< p(T) and f(T) 

In Fig||(A) the tree- level masses in eq. fl29| ) and m 2 (T) 
are shown for mP a eak (T = 0) = 550 MeV. m^(T) is not 
tachyonic and approaches to m 2 (T) in the symmetric 
phase. This confirms that our resummation procedure 



cures the problem of tachyons in Sec. II A. 

The solid line in Fig.||(B) shows the chiral condensate 
£(T) obtained by minimizing the effective potential with 



m 



peak 



(T = 0) = 550 MeV. f (T) decreases uniformly as T 
increases, which is a typical behavior of the chiral order 
parameter at finite T away from the chiral limit. 

As we approach the chiral limit (h — > or equivalcntly 
— > 0), £(T) develops multiple solutions for given T, 
which could be an indication of the first order transition. 
The critical value of the quark mass m" lt " below which 
the multiple solutions arise is 



m 



Tit./ Dhys. 



(mrv< hys f 



0.08, 



(38) 



where we have used Gell-Mann-Oakes-Renner relation 
26[ to related the pion mass with the quark mass, 
is the physical light-quark mass corresponding 



to mP^ 8 - 



140 MeV. The critical temperature for 



/??' 



crit. 



/mf ys - 



0.08 is T c 



170 MeV. The behavior 



of £(T) for m 7r (T = 0) = 30 MeV (just below the critical 
value m" 1 ' ) is also shown by the dashed line in FigJ|(B) 
for comparison. 



E. Chiral limit 

In Fig.|| the chiral condensates are shown for the chiral 
limit (m n = MeV) and for m, = 10 MeV. The phase 
transition looks like a first order in these cases. The 
existence of the multiple solutions of the gap equation for 
the 0(4) a model in the mean-field approach has been 
known for a long time ||27|j . Our analyses confirm this 
feature within the framework of OPT. 

However, this first order nature is likely to be an arti- 
fact of the mean-field approach as discussed in the second 
reference in |^7j : the higher loops of massless tt and al- 
most massless a are not negligible near the critical point, 



and they could easily change the order of the transition. 
In fact, the renormalization group analyses as well as the 
direct numerical simulation on the lattice indicate that 
the 0(4) a model has a second order phase transition 



We have also studied the free energy as a function of 
T near the chiral limit and found that it has a discon- 
tinuity near the critical point. This is another sign that 
the first-order nature is an artifact of the approximation. 
(Remember that, the free energy must be a continuous 
function of T irrespective of the order of the phase tran- 
sition.) 



F. Spectral function at finite T 

As one of the non-trivial applications of OPT, we show, 
in Fig.Q, the spectral functions of tt and a at finite T 
defined in ( p6|) . 

The figure shows that the spectral function of a, which 
does not show a clear resonance at T = 0, develops a 
sharp enhancement near 2ir threshold as T approaches 
T c . This is due to a combined effect of the partial restora- 
tion of chiral symmetry and the strong a — 2tt coupling. 
This is an typical example of the softening (or the precur- 
sor of the critical fluctuation) associated with the partial 
restoration of chiral symmetry [ p9| . The experimental 
relevance of this softening has been examined in the con- 
text of the low- mass diphoton production 0] . 

In the 7r-channel, a continuum develops by the scatter- 
ing with thermal pions in the heat-bath; tt + 7r thernial — > 
a. Because of this process, the pion acquires a width 
~ 50MeV at T = 145 MeV, while the peak position does 
not show appreciable modification. They are in accor- 
dance with the Nambu-Goldstone nature of the pion. 



IV. SUMMARY 

We have shown that the optimized perturbation theory 
(OPT) developed in J| naturally cures the problems of 
the naive loop-expansion at finite T, namely, the break- 
down of the naive perturbation at T > T c (due to the 
hard thermal loops) as well as at T < T c (due to the 
tachyonic poles). 

Furthermore, OPT has several advantages over other 
resummation methods proposed so far: 

First of all, the renormalization of the UV divergences, 
which is not a trivial issue in other methods, can be car- 
ried out systematically in the loop-expansion in OPT. 
This is because one can separate the self-consistent pro- 
cedure (Step 3 in Sec. II C) from the renormalization 
procedure (Step 2 in Sec. II C). 

Secondly, the Nambu-Goldstone (NG) theorem is ful- 
filled in any give order of the loop-expansion in OPT for 
arbitrary N in O(N) (f> 4 theory. This is because OPT 
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preserves the global symmetry of the effective potential 
in each order of the perturbation. 

There are many directions where OPT at finite T may 
be applied. The phase transition in 0(4) <p A theory as 
well as the dynamical critical phenomena near the crit- 
ical point are one of the most interesting problems to 
be examined further. The PMS condition for the effec- 
tive action will be suitable for this purpose. OPT may 
also have relevance to develop an improved perturbation 
theory for gauge theories in which the weak coupling ex- 
pansion is known to break down in high orders pOl . 
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(A) (B) (C) 

FIG. 1. Bubble and cactus diagrams. 






(A) (B) (C) 

FIG. 2. Diagrams which contain UV divergences as a result of the multiple insertion of (l/2)x<t> 2 - (A) corresponds to a single 
insertion with two external lines. (B) and (C) have no external lines with a single insertion and a double insertion, respectively. 
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FIG. 3. Spectral functions at T = in the n channel (A) and in the a channel (B) for m% eak = 550 MeV, 750 MeV and 1000 
MeV. 
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FIG. 5. (A) Masses in the tree-level moniT) and mo CT (T) shown with left vertical scale, and the mass parameter m 2 (T) with 
the right vertical scale. (B) f(T) for m v {T = 0) = 140MeV and 30 MeV with ml eak {T = 0) = 550 MeV. 
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FIG. 7. Spectral function in the tv channel (A) and in the a channel (B) for T = 50, 120, 145 MeV with m p / ah (T = 0) = 550 
MeV. 
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